Bidirectional transport on a dynamic lattice 
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Bidirectional variants ol stochastic many particle models for transport by molecular motors show 
a strong tendency to form macroscopic clusters on static lattices. Inspired by the fact that the 
microscopic tracks for molecular motors are dynamical, we study the influence of different types 
of lattice dynamics on stochastic bidirectional transport. We observe a transition toward efficient 
transport (corresponding to the dissolution of large clusters) controlled by the lattice dynamics. 



Microscopic models of stochastic transport like, e.g., 
the so-called asymmetric exclusion process (ASEP) or 
the zero range process (ZRP) have been extensively dis- 
cussed in the past two decades [l| . The interest is partly 
due to the fact that these models play a key role in devel- 
oping a general framework for statistical physics far from 
equilibrium . Moreover these models have been applied 
to many different transport problems, like road traffic Q. 
More recently variants of the ASEP have been used in or- 
der to describe intra-cellular transport phenomena driven 
by molecular motors [1, i-e., proteins that are able to 
transport cargos along the filament network of biological 
cells Q . In these models particles can also attach to and 
detach from the filament, leading to a finite path length 
of the molecular motors i, i, 0]. Remarkably, if open 
boundary conditions are applied, this extended model is 
able to predict the experimentally observed bulk localiza- 
tion of high and low density domains Q although many 
structural aspects of the filaments and motors have been 
neglected. 

In contrast to the success of these models describing 
the unidirectional collective motion of molecular motors 
in motility assays, it is still an issue to understand the rel- 
evant mechanisms involved in intracellular bidirectional 
transport. Several models of bidirectional transport have 
been suggested that consider positional exchange of par- 
ticles with different moving directions on the same track 
(see, e.g., [§] and [13 )• However, this family of mod- 
els cannot be used to describe the motion of oppositely 
moving molecular motors which cannot permeate each 
other on the same track. This scenario is consistent with 
recent findings for certain members of the kinesin and 
dyncin supcrfamily, moving in opposite directions on mi- 
crotubules (MTs, which arc polar filaments and a part of 
the cytoskeleton) and sharing the same binding site [ll| . 



Another difficulty arises from the fact that in real cells 
the volume surrounding the molecular tracks, i.e., the 
cytoplasm, is finite. Confinement introduces a kind of 
memory effect for the motors that have detached from the 
filament, i.e., they are more likely to attach again in the 
vicinity of their detachment position. As a consequence, 
domains of high particle densities on the track and in the 
diffusive surrounding reservoir co- localize [1, [T^. Both 
effects, i.e., inability of exchanges on the track and the 
explicit memory of the particle position, have been con- 
sidered recently for static one- and multi-lane systems. 
It has been shown [l^l that for generic model parame- 
ters the formation of macroscopic clusters and therefore 
a breakdown of the flow on the filament is observed. Re- 
markably the formation of stable macroscopic clusters 
does not depend on the particle density but rather on the 
number of particles in the system. Therefore the natu- 
ral question arises: What are the minimal prerequisites 
in order to maintain bidirectional stochastic transport of 
interacting particles in small volumes when they cannot 
cross each other on the same track? One suggestion has 
been made by Klumpp and Lipowsky [Tsl who considered 
that particles prefer to attach themselves in the neighbor- 
hood of motors of the same type, an effect that has been 
observed experimentally in vitro. For high densities and 
large differences in the binding affinity this leads to a 
spontaneous formation of unidirectional traffic. 

In the present paper, we propose a completely different 
type of mechanism which could lead to efficient bidirec- 
tional transport on a single track through consideration 
of the filament dynamics. This extension of the model has 
been inspired by the experimentally observed dynamics 
of the cytoskeleton. Indeed, the MTs on which molecular 
motors move are themselves highly dynamic IJ], due to 



" ebbinghaus@lusi . unl-sb . de] 



cecile. appert-rolland@th. u-psud.fr[ author 
dence should be addressed 
•l |l.santen@mx.uni-saarland.del 



to whom 



corrcspon- 



nucleation, polymerization and depolymerization, which 
occur on time scales similar to those involved in motor 
transport and are thus likely to interfere with the motor 
dynamics. Beyond the interest for intracellular traffic, 
we shall give here a first example where a dynamically 
driven jammed phase is hindered by the lattice dynamics. 

The model [T^ consists of two species of particles mov- 
ing on two parallel lanes (Fig. [T]) with periodic boundary 
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FIG. 1: Schematic representation of the particle dynamics 
(lattice dynamics is not included in this figure) . Arrows indi- 
cate possible moves with corresponding rates which are sym- 
metric for both particle species. We impose hard-core inter- 
action on the lower lane (filament), while the particles on the 
upper lane are non-interacting. Periodic boundary conditions 
are considered. 



conditions. Along the lower lane, which mimics the po- 
lar filament, the particles perform directed motion (rate 
p) in the direction determined by the particle's species. 
In the upper lane, particles diffuse freely (rate D) and 
do not interact (sites on the upper lane can be multiply 
occupied). Attachment to the lower lane happens at rate 
LUa- The detachment rate uJd is chosen to be smaller than 
the stepping rate p in order to capture the processivity of 
molecular motors. This setup is motivated by the quasi- 
one-dimensional geometry of axons although it largely 
(over-)simplifies the real structure. Please note that one 
can account for the geometry of the diffusive reservoir 
by tuning the attachment rate of the particles • Here 
we consider only one processive lane. Counterintuitively, 
this simplification reduces, for a given particle density, 
the tendency to form macroscopic clusters as shown in 

As a new feature we add some lattice dynamics for 
the lower lane, i.e., some sites are eliminated and recre- 
ated. The diffusive upper lane remains unchanged. On 
the lower lane, particles can only occupy a site if this 
binding site exists. The attachment moves (rate uja) are 
consequently rejected if the binding site has been elimi- 
nated. Additionally, a particle will automatically switch 
to the upper lane if it makes a forward step (rate p) onto 
an eliminated site or if the site which is occupied by the 
particle is eliminated. 

In the following, we consider different types of lat- 
tice dynamics. In particular: (1.) a simple realiza- 
tion of uncorrelated lattice dynamics. (2.) Dynamics 
which depends on occupation by particles. The results 
shown were obtained from Monte Carlo simulations over 
at least 10^ steps with a constant set of parameters for 
the particle dynamics: p = 1, uJd = 0.02, uja = 0.33 and 
D = 0.33. The particle density has been chosen high 
enough (p^^ = 0.5) in order to observe large clusters in 
the case of a static lattice even for small system sizes. 

First we consider some lattice dynamics independent 
of the configuration of particles: a site of the lower lane 
is randomly eliminated at rate kd and recreated at rate 
kp. This choice of the lattice dynamics can be seen as 
a minimal model, with a limited number of parameters. 
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FIG. 2: (color online). Distribution of cluster sizes in a system 
of size L = 1000 with simple lattice dynamics (scenario 1). 
Recreation of lattice sites occurs at fci, = 1. The black line 
con 
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FIG. 3: (color online). Flux of positive particles along the 
filament with simple lattice dynamics at fcp = 1 for different 
system sizes L and for the same density. If depolymerization 
is weak, the flux depends on the system size, but changes to a 
density-dependent state as large clusters disappear (see also 
Fig. [2} . The inset shows the same data divided by the flux in 
the smallest system (L = 100). 



Network dynamics induces an increase (resp. decrease) 
of the effective detachment rate Wd,cff (respectively, at- 
tachment rate Wo^eft), which take the form 
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As the depolymerization rate kd increases - and accord- 
ingly the fraction of "holes" in the filament kd/{kp + kd) 
- large clusters become less and less dominant until they 
disappear completely (Fig. [2]). This disappearance of 
large clusters is accompanied by a transition from a size 
dependent to a size independent state (Fig. [3]). For large 
systems, fiux vanishes below the transition, while it keeps 
a finite value above the transition. As a result, the tran- 
sition becomes sharper for an increasing system size. 

When the depolymerization rate kd is increased, first 
the flux along the filament (symmetric for both particle 
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FIG. 4: (color online). Space-time plots of the filament oc- 
cupation for a static (left) and a dynamic (right) filament in 
a system of size L — 1000 and density = 0.5. Red resp. 
green dots are particles with moving direction to the right 
resp. left, white spaces are unoccupied filament sites. The 
filament dynamics clearly induces a transition from a conden- 
sated to a homogenous phase. Note that considering a finite 
capacity of the diffusive lane would reinforce the jamming. 
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FIG. 5: (color online). Maximum flux j'(,(fcp, fed, max(fcp)) (red 
squares) and optimal depolymerization rate kd,n-iiix{kp) (black 
circles) at given polymerization rates kp in a system of size 
L — 1000. For comparison, the flux in a system with a static 
lattice is drawn as a solid red line. Note that the horizontal 
axis is in logarithmic scale and, therefore, both quantities 
seem to saturate at high polymerization rates. 



species) increases too. However, note that, if kd is too 
high, binding sites become increasingly sparse, and par- 
ticles have fewer segments on which they can contribute 
to the total flux in the system. Hence, the flux of each 
particle species along the filament disappears for kd ^ oo 
and kp constant. This behavior entails the existence of 
an optimal value for k^ (not seen on Fig. [3]) at which 
the flux is maximized. In this optimal regime with not 
too many holes in the filament, blocking situations with 
at least two particles of different species still occur fre- 
quently. As a consequence, the maximum fluxes that we 
obtained are about one third of the flux in a comparable 
single-species system [l^, which is nevertheless a great 
improvement compared to the static filament case. 

The degree of homogeneity can be estimated by com- 
parison with a product-state. Indeed, above the tran- 
sition, particles are well-dispersed over the whole sys- 
tem (Fig. [4]) and are much less correlated than in the 
static filament case. In a product-state approximation 
which neglects any correlations and if translational in- 
variance is assumed, the stationary state leads to [Hj 
'^a,off/o^(l - - pf) = Wd,offP^ where p denotes the 
particle density, i.e., the number of particles divided by 
the system size. The indices u and b refer to the un- 
bound state (upper lane) and bound state (lower lane) 
and +/— signs denote the particle species moving either 
in positive or negative direction. In combination with 
the conservation of particles = Ph^ + i the density 
of particles on the filament pb as well as the flux of one 
particle species along the filament = pp^{l — p^ ^ Pt) 
can be calculated. 

Below the transition, a big cluster is formed and 
the product-state approximation obviously fails. Above 
the transition, the product-state prediction still overesti- 
mates (not shown) the flux. This indicates that clusters 
made of a few particles still play an important role. It is 



only for strong lattice dynamics that the product-state 
solution regains validity. 

In Fig.[5l simulation results for the optimal depolymer- 
ization rate kd^uia.xikp) maximizing the flux are drawn as 
a function of the polymerization rate. The maximum 
flux jb{kp, fcci,max(fcp)) IS shown as a function of kp. One 
can see that the gain in transport capacity is obtained 
for a wide range of values of the polymerization rate kp - 
an observation which supports the general validity of our 
results. Furthermore, the fraction of eliminated filament 
sites kd/{kp -\- kd) in the optimal regime decreases with 
increasing polymerization rate kp as the optimal value 
kd,max saturates. This indicates that the optimal fiux 
depends on the time interval a site persists rather than 
on the delay after which it is recreated. 

In the second realization of lattice dynamics consid- 
ered, a site is eliminated with rate kd only if that site is 
occupied. This can be considered as a prototypical ex- 
ample of a coupling between MT dynamics and transport 
by molecular motors. On average, the fraction of existing 
sites on the lower lane becomes kp/ [kp + [p'^ + p^)fc(i] 
and Eq. ^ is modified accordingly. 

It turns out that the results are qualitatively the same 
as obtained for the first type of dynamics. The main 
difference is that the fiux in the homogenous phase is 
higher than in the first scenario. In this phase, when a 
motor is moving freely, it encounters less holes than in the 
case of scenario 1, as empty sites cannot be eliminated 
anymore. Proccssive runs along the filament are thus less 
frequently interrupted. By contrast, in the condensated 
phase, there is almost no difference to the first scenario 
as the transport capacity is limited by the large clusters, 
where all filament sites are occupied by particles anyway. 
The transition to a density-dependent efficient state is 
maintained if one generalizes scenario 2 by requiring that 
a higher number of particles has to accumulate in order 



4 



to put enough strain on the filament to break. 

Finally we mention briefly another kind of dynamics. 
Biopolymers, such as actin filaments or MTs, show a 
characteristic type of dynamics under certain conditions, 
termed tread- milling [6|, for which both ends of the fil- 
ament move with the same average velocity. We have 
checked that, for a fairly simplified model consisting of 
regularly spaced holes in the lower lane which propagate 
synchronously but stochastically through the system, the 
aforementioned transition toward efficient transport is 
still observed. Note that, since the holes move only in 
one direction, the two species of particles are affected 
differently. The flux of the two species therefore depends 
on the moving direction along the filament. For both 
species we remark a considerable increase of the flux and 
reach a maximum current comparable to the one in the 
previous scenarios. 

To summarize, we have presented a model for bidi- 
rectional transport on a one-dimensional dynamic lattice 
coupled to a confined diffusion reservoir. While on a 
static lattice persistent clusters form which inhibit effi- 
cient transport |12| . we have found the counterintuitive 



effect that the transient suppression of filament sites dra- 
matically enhances the transport capacity. Indeed, the 
inhibition of large clusters leads to a transition toward 
a homogenous state characterized by efficient transport. 
This is a new mechanism in the phenomenology of dy- 
namic phase transitions. This transition separates a size- 
dependent (jammed) state from a density dependent one. 
It is robust in the sense that it is obtained for quite dif- 
ferent types of lattice dynamics. The actual transport 
capacity of the system rather depends on the optimal 
lifetime of a binding site than on the details of the fila- 
ment dynamics. The lifetime has to be short enough in 
order to avoid jam formation and long enough in order 
to direct the transport. 

Some insight could be gained from an analysis simi- 
lar to the one found in [l7| for the symmetric motion 
of non-interacting particles in fluctuating energy land- 
scapes. However, note that here the effective "potential" 
landscape emerges spontaneously from the particle jam- 
ming. 

Finally we would also like to discuss the relevance of 
our model results for axonal transport. In view of our re- 
sults for bidirectional transport on static lattices, which 



generically leads to jamming, it is rather surprising, but 
of course necessary, that transport in these real one- 
dimensional axonal structures would be at all stable and 
efficient. This is all the more surprising given that the 
model overestimates the diffusivity of the detached par- 
ticles, which should be reduced in real systems due to 
the size and interactions of the vesicles in the cytoplasm. 
These results give strong evidence that additional effects 
must come into play in order to stabilize motor driven 
transport in axons and reduce the tendency to form large 
particle clusters. The mechanism suggested in this work 
is based on the fact that the filament dynamics limits 
the size of particle clusters. A more detailed modeling 
of motor driven transport in axons is difficult since the 
experimental findings concerning the motor interactions 
and the MT-dynamics are so far incomplete and subject 
to interpretation: while it is well established that the 
plus-ends of MTs undergo polymerization events toward 
the synapse, the pathlength of the growing plus-ends as 
well as the dynamics of the minus ends is not known yet. 
Besides the robustness of our results with respect to the 
lattice dynamics details, the importance of the lattice dy- 
namics is supported by experimental results on transport 
in axons, where the strong correlation between MT dy- 
namics and vesicle transport has been demonstrated [l^ ■ 
In this context, we have observed that lane formation, 
which could be invoked as an alternative mechanism for 
bidirectional transport is impeded by the dynam- 
ics of MTs. Thus it seems that it is necessary to take 
into account the lattice dynamics in order to gain un- 
derstanding in axonal transport. Although our picture 
of bidirectional transport is qualitatively consistent with 
the experimental findings it is still an issue to develop a 
more quantitive description of real axonal transport. A 
fully validated model could then be used in order to inves- 
tigate mechanisms that are underl y in g clinically relevant 
transport defects in nerve cells []J,ll9l. 
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